Nuclear Magnetic Orderings and Frustration in Bcc 3 He in High Magnetic Fields 
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We studied phase transitions and thermodynamic properties of the field-induced antiferromagnetic 
nuclear-spin ordered phase in bcc solid 3 He by measuring the melting pressure at temperatures down 
to 0.4 mK in magnetic fields up to 15 T. The transition temperature from the paramagnetic phase 
is reentrant with increasing field with a maximum at 10 T. This indicates that the system is highly 
frustrated by the competing multiple-spin exchanges and excludes other mechanisms related to the 
zero-point vacancies. The upper critical field was estimated as 19.7 T, which reveals non-negligible 
contributions from higher order exchanges beyond six spins. A considerable softening of the spin 
waves in the ordered phase in low fields also suggests the strong frustration. 
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Solid 3 He with bcc structure is an ideal three di- 
mensional antiferromagnet with intrinsic frustration due 
to the competing multiple-spin exchanges (MSEs) 0. 
Nearly perfect bcc crystals can easily be grown when 
they coexist with the liquid phase. It has a two dimen- 
sional (2D) counterpart, monolayer solid 3 He adsorbed 
on graphite, where an exotic ground state so called the 
gapless quantum spin liquid is realized 0. Thus we can 
study the frustration due to the MSEs and its dimen- 
sionality effects in great details through investigations of 
nuclear magnetism of bcc and 2D 3 He. 

The exchange interactions among the nuclear spins 
(S — 1/2) in solid 3 He are associated with the purely 
isotropic direct atom-atom exchanges. They completely 
dominate the thermodynamic properties below 100 mK, 
which makes this material an exceptionally pure mag- 
netic system. As Thouless first conjectured, the frus- 
tration in this system originates from the competition be- 
tween antiferromagnetic (AFM) and ferromagnetic (FM) 
MSEs. He showed that even numbers of exchanging 
atoms introduce AFM interactions while odd numbers 
do FM ones. The importance of the four-spin ring ex- 
change has recently been recognized in highly correlated 
electronic systems as well 0. 

It is a remarkable and unique feature of solid 3 He that 
even absolute values of the MSE frequencies (Jp) are 
potentially calculable from the first principles. Jp up 
to six-spin exchanges have been eventually calculated by 
the path integral Monte Carlo (PIMC) technique [j| for 
bcc 3 He at a molar volume (24.12 cm 3 /mol) close to the 
melting one. They are consistent with most of the so far 
existing experiments semi-quantitatively (within 30%) |(| 
7]. It is, however, generally difficult to test them more 
rigorously because of the many parameters involved. In 
addition, theoretical approximations applied to the MSE 
Hamiltonian sometimes limit the accuracy of the tests. 
Therefore, there still remains a fundamental question for 
the MSE model; Where can we truncate the series of the 



MSE interactions? 

The competing interactions in bcc 3 He might be ex- 
plained by quite different hypotheses from the MSE 
model. The most intriguing one is so called the zero- 
point vacancy (ZPV) model 8]. If the system contains a 
finite amount of ZPVs, FM interactions will be induced 
among nuclear spins surrounding the ZPVs in order to re- 
duce their kinetic energies. Such FM interactions would 
compete with the ordinary AFM two-spin exchange in- 
teraction. This model can be quantitatively tested by 
exploring the high-field magnetic phase diagram at the 
melting density where the ZPV concentration is maxi- 
mized. This test is of up-to-date interest since ZPV could 
be responsible for the superfluid responses recently ob- 
served in solid 4 He |9|. It is also important in the light 
of determining the zero-temperature upper critical field 
B C 2{Q), beyond which the spins are fully polarized along 
the external field, for this field is exactly calculable even 
with the mean field theory . 

In this Letter, we present results of melting pressure 
measurements of bcc 3 He down to temperatures well be- 
low the AFM ordering temperature (T c ) and in high mag- 
netic fields up to 15 T, a factor of two higher field than in 
the previous experiment |l(J. The melting pressure (P) 
is a useful probe to study thermodynamic properties and 
phase transitions in 3 He. The temperature variation of 
P is determined by the Clausius-Clapeyron relation: 



dP/dT = (Si - S)/(Vi - V), 



(1) 



where S (Si) and V (Vi) are the molar entropy and vol- 
ume of the solid (liquid) phase, respectively. Note that, 
below 5 mK, Si is less than 10% of S and precisely known 
from the specific heat measurement 111. The volume 
difference (Vi — V) — 1.314 cm 3 /mol [l2| is expected to 
be temperature and field independent. Thus the melt- 
ing pressure change is determined predominantly by the 
entropy of the solid phase. Details of the experimental 
setup have already been described elsewhere 0. We 
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FIG. 1: (a) Entropy, (b) specific heat and (c) dP/d(T 4 ) of bcc 
3 He deduced from the melting pressure (P) data at B = 2.0, 
7.9 and 14.5 T. The vertical arrows indicate the HFP-PP tran- 
sitions. The solid lines are guides to the eye. The open circles 
are from the previous direct specific-heat measurements [l4|. 



took all the data under complete thermal equilibrium. 

Figure 1 (a) (b) show the entropy and the specific heat 
(C) of bcc 3 He deduced from the measured melting pres- 
sure through Eq. at three selected fields. The HFP- 
PP transitions are identified as kinks in the S vs. T plots 
and as peaks in the C vs. T plots. These features in- 
dicate that the transitions are not first order but con- 
tinuous. This is consistent with the direct specific heat 
measurements below 8 T [3 ^| , and with the theoret- 
ical proposal that the spin structure of HFP is canted 
normal antiferromagnetic (CNAF) with two sublattices 
and cubic symmetry 0] (see Fig. 2). T c can also be 
determined as a temperature where dP/d(T 4 ) has a min- 
imum in the plot as a function of T 4 (Fig. lc). The T 4 
dependence of P(T) is expected from the AFM spin wave 
theory in the ordered phase. The above three different 
determinations of T c agree each other within 5% except 
for 2 and 6 T (13 and 16%, respectively). 

Figure 2 shows the magnetic phase diagram of bcc 3 He 
at 24.21 cm 3 /mol. The data points denote phase bound- 
aries between two of the three magnetic phases (HFP, 
LFP and PP) determine d by the present and previous ex- 
periments [E E3, EE USUI Here the low-field phase 
(LFP) is another AFM ordered phase with four sublat- 
tices (the U2D2 phase ^^) which exists below the lower 
critical field [B c \ = 0.45 T). Note that we made small 
volume corrections for the original data so as to be con- 
sistent with the result for the fixed V at 24.21 cm 3 /mol 
in Figs. 2, 3(b) and 4. This was necessary because V 
and hence Jp increase slightly with increasing B along 
the melting curve. The largest correction is 4.6% at 
B = 14.5 T. 

Obviously the HFP-PP boundary has a reentrant 
shape with a maximum T c (= 3.4 mK) near B = 10 T. 
This is the first experimental evidence that T c turns 
around to decrease above a certain field at the melt- 
ing density as was observed at a much higher density 
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FIG. 2: Nuclear magnetic phase diagram of bcc 3 He at 24.21 
cm 3 /mol. The data points are phase transitions between two 
of the three phases (HPF, LFP and PP); this work (•), Ref. 
E3 (o), Ref. (□), Ref. [T^ (A) and Ref. (+). Also 
shown are the HFP-PP transitions at 22.47, 22.69 and 23.06 
cm /mol from Ref. HE3- 

The solid and dash-dotted lines 
are HFP-PP boundaries scaled from that determined for 22.69 
cm 3 /mol with F = 18.5 and 21.3, respectively. For the other 
lines see the text. 

[6J. From this, we conclude that the HFP is an AFM 
ordered phase and not any phase with high polariza- 
tion associated with the ZPV. Within the ZPV model 
0, upon increasing magnetic field, the vacancy band 
width increases due to the increasing polarization and 
the vacancy creation energy decreases due to the decreas- 
ing melting pressure. Hence the vacancy concentration 
n should monotonously increase, which never yields the 
negative slope for the phase boundary. Generally, the 
Griineisen constant for n should be quite different from 
that for J P , i.e., T(J P ) = d\nJ P /d\nV « 18 7]. Nev- 
ertheless, the measured HFP-PP boundary has a similar 
volume dependence to L(Jp) as the data points follow 
approximately the solid line estimated from the bound- 
ary at 22.69 cm 3 /mol @ with L = 18.5 (see Fig. 2). 

The dotted line in Fig. 2 is a CNAF-PP boundary cal- 
culated by the Green function method |20J with the fol- 
lowing MSE parameters: 

J 1N = 0.48, J 2N = 0.067, Ti = 0.20, K P = 0.28, 
K F = 0.028, K A = 0.006, K B = 0.0005, K L = 0.011, 
K s = 0.0020, Si = 0.037, S 2 = 0.023, (2) 

given by the PIMC calculation 5]. Here, all numbers are 
in mK, and we followed Ref. [2l| for the notations of Jp; 
two-spin (Jitv, J2iv), three-spin (Ti), four-spin [Kp, Kp, 
Ka, Kp, Kp, Ks) and six-spin (Si, S2) exchanges. Note 
also that all the original MSE parameters in Ref. |5j have 
been increased only by 4% in Eq. @ so that they give 
B C 2(0) — 19.7 T, the estimation discussed later. The 
agreement between the Green function calculation and 
the present data is good. The dashed line is the HFP- 
PP boundary calculated by the mean field approximation 
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FIG. 3: (a) Log-log plots for the zero-temperature upper and 
lower critical fields, B&iO) and B c i(0), as functions of V; 
Ref. H E3 (•), Ref. [H 13 (A), Ref. ^ (■), Ref. E3] 
(♦), estimations from the HFP-PP boundary (o) and Eq. Q 
(0)- The dashed lines are MFA calculations with the MSE 
parameters of Eq. (|5J and their volume dependencies given 
in Ref. [HEil. (bVMagnetization of bcc 3 He at T = 0; this 
work (•) and Ref. ^ (A). The solid line represents Eq. @ 
and the thin line is the linear extrapolation of the highest 
three data points. The dashed line is the MFA calculation 
with Eq. The dotted line is a guide to the eye giving 

B c2 (0) = 19.7 T. 



(MFA). It gives too high T c compared to the experiment 
particularly at high fields indicating the importance of 
thermal fluctuations @. 

An accurate estimation for B c2 (0) at the melting den- 
sity is available by extrapolating the volume dependence 
of this quantity which is known more reliably at higher 
densities. As is shown in Fig. 3(a), all the existing data 
[E El HE H3 follow the power law dependence with 
r = 19.9 ±0.2 very well, which gives 19.7 ±0.4 T at 24.21 
cm /mol. This Griineisen constant is considerably larger 
than those for any other physical quantities measured 
at lower fields (« 18). The WKB calculations [HUH 
predict the larger T for the smaller Jp in the density 
range for bcc 3 He, i.e., T(J 1N ) = 18.2, T(Ti) = 16.2, 
T(Kp) = 15.1, T(K F ) = 18.1, and F(Si) = 18.9 which 
explain very well the exceptionallysmall experimental T 
value (= 15.2 ± 0.2) for B c i(0) B E3 0n the other 
hand, the same WKB calculations predict T = 18.1 for 
B c2 (0) which is too small compared to the experiment 
(see Fig. 3(a)). It is important to recognize that higher 
order exchanges contribute more effectively to magnetic 
properties at higher fields particularly to B c2 (0) [M l2l| . 
The discrepancy can not be explained even by infinitely 
large T values for Si and S 2 in Eq. ©. We thus conclude 
that the convergence of the series of MSE interactions is 
somewhat slower than the current assumption and that 
higher order exchanges than six spins should be taken 
into account. If so, at fields above 15 T, the B C 2(T) 
line would deviate upward from the solid line in Fig. 2, 
which assumes the constant T value (= 18.5), in order 
to smoothly connect to 19.7 T at T = 0. Note that the 
ZPV model gives a lower T value than 18 because of the 



increasing FM contribution at larger volumes, which is 
inconsistent with the present result. 

It is intriguing to note that T c approaches zero as B — » 
(inset of Fig. 2). In other words, if the LFP did not 
exist hypothetically, the system would have no long range 
order at B = down to zero temperature due to the 
strong frustration. Actually, this is realized in 2D 3 He 
|2|. Thus a similar high- field ordered phase may exist in 
2D as well 26] . Another point is that even the highest T c 
(= 3.4 mK) is five times lower than the exchange energy 
represented by [iB c2 (0) (— 16 mK), where /i is the nuclear 
magnetic moment of 3 He. This large suppression of T c is 
a further evidence for the strong frustration inherent in 
this system. 

We extrapolated the melting pressure data below 0.6 
T c to T = by fitting them to 



P(T,B) = P(0,B)-cT 4 



(3) 



The zero temperature magnetization M was deduced 
from the fitted P(Q, B) using the mag netic Clausius- 
Clapeyron eq. as was done in Ref. [27j . The reduced 
magnetization, m = M/M sat , obtained in this way is 
shown in Fig. 3(b). Here M sat is the saturation magne- 
tization. Our data can be fitted well to the thick solid 
line expressed by 

B = -(10.1 ± 1.7)m + (25.4 ± 5.2)m 3 + (6.2 ± 3.9)m 5 (4) 

with an rms deviation of 1.5%. This functional form is 
an MFA expression for the magnetization curve based on 
the MSE Hamiltonian considering up to the six-spin ex- 
changes [2l] . A simple extrapolation of Eq. Q to m — 
1 gives B C 2(0) = 21.5 ± 0.4 T, which agrees well with 
the similar extrapolation in Ref. (= 21.7 ± 1 T) but 
is higher than the above mentioned estimation (= 19.7 
T) by 9%. This is an additional indirect evidence for 
the contributions from the currently ignored higher or- 
der exchanges. The magnetization seems to have a weak 
positive curvature above 15 T in order to be consistent 
with B c2 (0) = 19.7 T H|. One should be careful about 
the accuracy of MFA, which ignores fluctuations, even 
near B c2 (0) since it gives too high B c2 at finite temper- 
atures. However, a better approximation taking account 
of fluctuations such as the coupled cluster approximation 
(CCA) gives essentially a similar magnetization curve to 
the MFA one [25|. 

From the fitting parameter c in Eq. J3J), we deduced the 
angle averaged spin-wave velocity (v) using the relation: 



Nir 2 k% 



V 



90ft 3 {Vi - V) v 3 ' 



(5) 



The number of spin-wave modes (N) is assumed as unity, 
since thermal energies at T < 0.6T C are much lower than 
the Zeeman energies in the field range studied here. In 
Fig. 4 we plot the deduced v as a function of B. The 
field dependence is qualitatively similar to that of T C (B), 
because T c is a measure of the spin-wave stiffness. The 



4 



24.21 cm 3 /mol 




o 1 1 ■ . , . i . . " , i . y-°. , i , ■, , , -i , , , , i 

5 10 15 20 25 

B (T) 

FIG. 4: Spin wave velocities in the HFP and LFP; this work 
(•), Ref. tOt (V), Ref. |2?| (A), Ref. (+), Ref. [U (o). 
The solid line is a cubic fitting of the present data. The dotted 
line is an SWA calculation with Eq. for the CNAF phase. 
The same calculation for the U2D2 phase is also shown (x). 



smooth extrapolation to lower fields agrees well with the 
previous data below 1 T [Hl27ll29j. 

Although our linear spin- wave approximation (SWA) 
calculation for v with the MSE parameters in Eq. 
(the dotted line in Fig. 4) agrees with the experiment 
at the highest fields, it gives too large values in lower 
fields, e.g., four times larger at B c \. In order to obtain 
such small v values at low fields, quite large J2N and S\ 
even comparable to Jin have to be assumed in the cal- 
culation, which is unreasonable. Thus the discrepancy 
indicates that the spin-wave excitation is ill defined in 



the HFP in low fields due to the frustration and the re- 
sultant large fluctuations. The fluctuations are expected 
to be larger in the CNAF phase than in the U2D2 phase 
because of the non-collinear spin structure in the former 
phase. Eventually, the same SWA calculation for the 
U2D2 phase yields v = 9.0 cm/s which is in reasonable 
agreement with the experimental one (= 7.7-7.8 cm/s) 
|29ll30| . To test this explanation, it is worthwhile to per- 
form the next-order spin wave calculations for both the 
ordered phases. 

In conclusion, we studied the nuclear AFM orderings 
of the melting bcc 3 He in high fields. The observed reen- 
trant nature of the HFP-PP boundary strongly support 
the MSE hypothesis excluding other scenarios related to 
the ZPV. The zero-temperature upper critical field has 
been determined as 19.7 T reliably from its volume de- 
pendence. The higher order exchanges than six spins 
can not be neglected to account for the present high field 
data. We also found the considerable softening of the 
spin waves in the HFP at low fields, which is indicative 
of strong frustration due to the competing MSEs. A fu- 
ture extension of this work to higher fields near 20 T and 
a new PIMC calculation of Jp beyond six spins at differ- 
ent densities are desirable to verify these conclusions. 
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